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were calculated using the full-potential LMTO method. The metal vacancies change the density 
of states near the Fermi level and this effect is quite different for III-V group transition metal 
diborides. Contradictory data on superconductivity in diborides may be supposed to be connected 
with nonstoichiometry of samples. Vacancy formation energies are estimated and analyzed. 
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A recent discovery of superconductivity (SC) in MgB2 
(T c w 39 Kp and creation of promising materials based 
thereon (see reviews' 2-5 !) have attracted a great deal 
of interest in hexagonal transition metal (M) diborides 
isostructural with MgB 2 ' 2-5 l. Earlier investigations ob- 
served no superconductivity in d-metal diborides (Ti, Zx, 
Hf, V, Ta, Cr, Mo, Nb) with temperatures above 0.6 KB. 
Experimental and theoretical investigations' 7-17 ! showed 
that the supercoductivity with T c « 40K is unlikely in all 
undoped diborides except MgB 2 , where the high value of 
T c is explained by strong electron-phonon coupling and 
high phonon frequencies and is closely connected with 
the existence of hole doped boron a bands. 

Recently, rather high T c have been reported for ZrB 2 
(5.5K)LI, TaB 2 (9.5K)i and_NbB2 (5.2K)H. On the other 
hand, Rosner et al showed!!] that-.SC in TaB 2 was absent 
down to 1.5 K, and according tooB the SC transition for 
NbB 2 is .not observed to T ~ 2K, in conformity with ear- 
lier dataQ. The weak electron-phonon coupling estimated 
to be inconsistent with superconductivity in ZrBji-3 leads 
to T c ~ 0.1 KpHji TaB 2 E3 and is responsible only for T c 
~ 3K in NbB 2 t 2 J. So, both experimental and theoretical 
studies performed to characterize and understand super- 
conductivity in MB 2 yielded contradictory results. It is 
important also to note that in all the studies' 7-17 ' the 
composition of transition metal diborides was considered 
to be strictly stoichiometric (B/M — 2). As is known, lat- 
tice vacancies are typical defects, and their presence may 
change considerably the properties of non-stoichiometric 
materials in the homogeneity region. For instance, one 
of the most familiar classes of non-stoichiometric com- 
pounds with exclusively wide homogeneity regions (to 30 
~ 55 at %) includes cubic (Bl) III-V group transition 
metal carbides, nitrides and oxides extensively investi- 
gated in numerous works, seaia. As distinct from the 
above Bl - carbides, nitrides and oxides, III-V group 
metal diborides have very narrow homogeneity regions 
under equilibrium conditionsEj. Therefore the role of 
non-stoichiometry is usually ignored when their proper- 
ties are examined. No theoretical studies describing the 
effect of lattice vacancies on the electronic structure of 
metal diborides are familiar to us either. 



Quite recently, Yamamoto et alE2l have carried out a high- 
pressure synthesis of a series of metal-deficient samples 
of the AlB 2 -type Nbi_ K B 2 , Tai_^B 2 (0 < x < 0.48) and 
Mo 2 B5-type borides with the compositions of MoB 2 and 
Moo.84B 2 .o- It was also found that for Nb 1 _ a; B 2 the su- 
perconductivity appears above x ~ 0.04. As vacancy con- 
centration increases, T c grows with the maximum value 
above 9K at x ~ 0.24. 

In this report, we are focusing on the effects of metal 
lattice non-stoichiometry on the electronic properties of 
metal diborides. For this purpose, the band structure 
of hexagonal NbB 2 , ZrB 2 and YB 2 containing 25 % M- 
vacancies was investigated theoretically. The choice of 
the phases allowed us to discuss the role of metallic va- 
cancies in modifying the band structure of isostructural 
MB 2 phases for all possible types of band filling of these 
materials. 

The above diborides have a hexagonal crystal structure 
(AlB 2 -type, space group P6/mmm) composed of layers of 
trigonal prisms of M atoms in the center of boron atoms, 
which form planar graphite-like networks. The metal- 
deficient phases Mo.75B 2 were simulated by 12-atom su- 
percells (2x2x1). Band structure calculations were car- 
ried out using a scalar relativistic full-potential linear 
muffin-tin method (FLMTO) with the generalized gradi- 
ent approximation (GGA) for correlation and exchange 
effects! 2 -! Latticg-,parameters of NbB 2 , ZrB 2 and YB« 
were taken fromE3. As follows from the experimented, 
the lattice constants for the metal-deficient Nbo.75B 2 (a 
= 3.098 A, c/a = 1.072) differ from those for the stoichio- 
metric NbB 2 by less than ~ 0.3 %. Therefore the lattice 
constants for the hypothetical Zro.7sB 2 and Yo.7sB 2 mere 
assumed equal to the a, c values of complete phasesE3. 
Calculated total densities of states (TDOS) and site- 
projected 1-decomposed DOS (LDOS) of NbB 2 , ZrB 2 and 
YB 2 and metal-deficient Nbo.7sB 2 , Zr .75B 2 and Y .75B 2 
are shown in Figs. 1 and 2, respectively. The band struc- 
tures of AlB 2 -like 4d-metal diborides are represented by 
bonding states formed by B2s and mixed B2p-Md states 
(A and B, Fig. 1) separated from the antibonding states 
(C) by a pseudogap. Depending on the M sublattice 
(the number of valence electrons n e ) , three types of band 
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filling of diborides are possible. For ZrB2 (n e =3.33 
e/atom), the Fermi level (Ep) is located in the pseu- 
dogap (TDOS minimum). This corresponds to the max- 
imum of chemical stability when the bonding states are 
completely occupied and antibonding states are vacant, 
see aW 23-25 l. When going to NbB 2 (n e =3.66 e/atom), 
the antibonding states become partially occupied and the 
TDOS at the Fermi level (N(E F )) increases, Table 1. On 
the contrary, for YB2 (n e =3.0 e/atom) some bonding d- 
p states are partially unfilled. This determines a decrease 
in the cohesive properties of NbB 2 and YB 2 as compared 
with ZrB 2 . Thisj-simple picture is consistent with the ex- 
perimental dataE3 and the first-principle calculations of 
the formation energies for these diborides, seel 23 ~ 25 l and 
below. 

Let us consider the main effects of M vacancies on the 
band structure of these diborides. The modification of 
electronic properties when going from MB 2 to MB0.75 
is generally determined by (i) a decrease in n e and (ii) 
changes of electronic states of atoms surrounding M va- 
cancies. According to the results obtained, these changes 
are quite different for III-V group metal diborides. Nb 
vacancies in NbB 2 result in the appearance of a new DOS 
peak and in a variation of DOS distribution in the region 
of the Fermi level (Figs. 2-4). The vacancy s-like states 
are below Ep. The modification of valence states near the 
Nb vacancy is visually demonstrated in Fig. 5. It is seen 
that no new Nb-Nb bonds going through the vacancy are 
formed. Deformation of Nb charge density contours takes 
place reflecting a growth of the electronic density along 
the Nb-Nb bond lines in the vicinity of the defect. For 
Nbo.75B 2 , N(E_p) decreases insignificantly (by ~ 1.9 %) 
relative to NbB 2 . 

In comparison with NbB 2 , for ZrB 2 the mentioned 
changes are much more pronounced, Fig. 2. The Zr va- 
cancies induce a new sharp DOS peak located in the pseu- 
dogap region. As a result, for metal-deficient Zro.75B 2 , 
N(Ej?) increases abruptly (from 0.30 to 1.22 states/(eV 
cell)). Quite different trends are observed in the band 
structure of Yo.7sB 2 : Ep is located in the local TDOS 
minimum, and N(Ef) reduces more than twice from 
0.900 to 0.404 states/(eV cell). The changes of LDOS 
contributions near the Fermi level from M and B states 
for III-V group metal diborides in the presence of M 
vacancies turned out to be different too, Table I. For 
Nbo.75B 2 , N(Ej?) diminished owing to a decrease in the 
contribution from Nb4d states (N Nb (Ep)), whereas the 
contribution from boron 2p states (N B (E F )) increases. 
On the contrary, a sharp growth of N(E^) for Zro.7sB 2 is 
due to growth of both N Zr (E F ) and N B (E F ), whereas a 
decrease in N(E^) for Yo.7sB 2 is caused by simultaneous 
reduction of N r (E F ) and N B (E jF ). 

Obviously, these changes in the DOS will affect the prop- 
erties of non-stoichiometric diborides, which depend, in 
particular, on the near-Fermi level electron density. Our 
estimations of the electronic specific heat coefficient (7, 
in the free electron approximation 7 =(7r 2 /3)N(E F )k 2 B) 
for complete MB 2 and metal-deficient M .75B 2 phases 



are listed in Table 2 as an example. It is seen that if 
for MB 2 7 changes in the order NbB 2 > YB 2 > ZrB 2 , 
for non-stoichiometric borides this sequence is as follows: 
Zr .7 5 B 2 > Nbo. 75 B 2 > Y .7 5 B 2 . The effect of M vacan- 
cies on the value of N(E^) is most noticeable for ZrB 2 . 
Thus, one can suppose that the value of T c ~ 5.5 K ob- 
served by Gasparov et apJor ZrB 2 may be due not only 
to the presence of ZrBi 2 E3 but also to nonstoichiometry 
in the Zr sublattice. 

For Nbo.75B 2 , no considerable changes are found near Ep. 
Obviously, the important role in the observed increase of 
T c for Nbi-^B^ will belong to the changes in phonon 
frequencies with lowering of crystal stability. On the 
contrary, a decrease in N(E^) for the non-stoichiometric 
Yi_2;B 2 does not allow us to expect superconductivity in 
the homogeneity region of yttrium diboride. 
Finally we discuss the variation in chemical stability and 
cohesive properties of non-stoichiometric diborides from 
the numerical estimations of the cohesive energy (E co h) 
and the heat of formation (AH) of MB 2 and Mo.7sB 2 
phases. The cohesive energies of MB 2 and Mo.7sB 2 were 
calculated as: 



coh 
pMo.75 -B: 
^coh 



nMB, 



^tot 

Tp Afo. 75 B-2 

- &tot 



{E%+2E B } 
■ {0.75E™ + 2E*}, 



(1) 



where E^ t , E B t are the total energies of free M and 
boron atoms, and E^ 82 , E^' 75S2 are the total energies 
(per formula units) of MB 2 and Mo.75B 2 , respectively. 
The heats of formation are defined as: 



ah mb 2 



MB 2 



= (E& + 2E? ot ) - E, 
(0.75E^ t + 2E t B ot )-E^r 5B2 , 



(2) 



where Ef^ t , Ef ot are the total energies of elemental 
metals and a-boron obtained from FLMTO calculations. 
In conformity with experimentslia, the MB 2 phases have 
positive values of formation heat. The value of AH MS2 
decreases (AH(ZrB 2 ) > AH(NbB 2 ) > AH(YB 2 )) being 
maximum for ZrB 2 , Table II. It should be noted that 
the cohesive energies (characterizing the atomic decom- 
position of MB 2 ) exhibit a different sequence: E^ B2 > 

^coh 2 > ^c.oh ■ The appearance of a vacancy sharply de- 
creases the crystal stability (AH, Table II). The energy 
of vacancy formation in the M sublattice was estimated 
as: 



E 



vf 



77 



MBa 



E Mo,,B 2 _ 0.25^, 



(3) 



The results obtained indicate (Table II) that the 
energy of formation of M vacancies in NbB 2 is lower 
than that in ZrB 2 and YB 2 . Thus, the M-sublattice 
nonstoiciometry is more difficult to achieve in ZrB 2 , 
which is the most stable diboride among those considered 
above. 

In conclusion, we have presented the first results of band 
structure calculations for metal-deficient Nb, Zr and Y 
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diborides performed by the FLMTO method. It was 
established that vacancies are more likely to appear in 
diborides of transition metals of III and V groups. The 
effect of M vacancies on the cohesion and electronic 
properties of III-V group transition metal diborides was 
found to be determined by the band filling. The changes 
in N(Ej?) caused by M vacancies are different: N(Ef) 



may decrease (YB2), increase (ZrB2) or stay almost 
constant (NbB 2 ). 
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TABLE I: Total (N(E F ) tot ) and partial (N(E F ) ; ) density of 
states at the Fermi level (states/eV/f.u.) and electronic spe- 
cific heat coefficients (7, mJ*mol _1 *K -2 ) for complete and 
metai-deficient Nb, Zr and Y diborides. 



Phase 


N(E F )' ot 


N(E F ) Ms 


N(E F ) AJp 


N(E F ) AJd 


N(E F ) A1/ 


N(E F )" S 


N(E F ) Bp 




7 Bil 


NbB 2 


f.Of 


0.00 


0.01 


0.65 


0.04 


0.01 


0.13 


2.39 


2.42 


Nb .75B 2 


0.99 


0.01 


0.02 


0.54 


0.02 


0.01 


0.15 


2.34 




ZrB 2 


0.30 


0.00 


0.00 


0.17 


0.01 


0.00 


0.04 


0.71 


0.67 


Zro.7sB2 


1.22 


0.01 


0.03 


0.43 


0.02 


0.02 


0.33 


2.87 




YB 2 


0.90 


0.01 


0.02 


0.36 


0.01 


0.00 


0.14 


2.12 


2.03 


Y0.75B2 


0.41 


0.01 


0.01 


0.12 


0.01 


0.00 


0.08 


0.96 





TABLE II: Heat of formation AH, cohesion energy E co h and 
vacancy formation energy of vacancy formation E„/ (Ry/f.u.) 
for stoichiometric and metal-deficient Nb, Zr and Y diborides. 



System 


E MB 2 
coh 


AH"'" 2 


T? M Q. 75-^2 
coh 




E„/ 


Nb - B 


1.82 


0.27 


1.57 


0.19 


0.08 


Zr - B 


1.72 


0.357 


1.46 


0.22 


0.13 


Y- B 


1.42 


0.26 


1.23 


0.16 


0.10 



■5 




FIG. 1: The total density of states of NbB 2 , ZrB 2 and YB 2 . 
The energies are relative to the Fermi level. 
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FIG. 2: The electronic density of states of Nbo.75B2.The en- 
ergies are relative to the Fermi level. 
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FIG. 3: The electronic density of states of Zro.7sB2. The 
energies are relative to the Fermi level. 




FIG. 4: The electronic density of states of Y0.75B2. The 
energies are relative to the Fermi level. 
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FIG. 5: Charge density on a hexagonal Nb plane in the vicin- 
ity of metal vacancy (V) in Nbo.7sB2. 



